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ABSTRACT 

The broad MgllA2800 doublet has been much studied in connection with its 
potentially important role as a virial estimator of black hole mass in high red- 
shift quasars. An important task is therefore identification of any line components 
likely related to broadening by non-virial motions. H igh s/n median com posite 



spectra (binned in the "4D eigenvector 1" context of ISulentic et al.l 120071 1 were 
constructed for the brightest 680 SDSS DR7 quasars in the 0.4 < z < 0.75 
range where both MgllA2800 and H/3 are recorded in the same spectra. Com- 
posite spectra representing 90% of the quasars confirm previous findings that 
FWHM(MgllA2800) is about 20% narrower than FWHM(H/3). The situation is 
clearly different for the most extreme (Population A) sources which are the high- 
est Eddington radiators in the sample. In the median spectra of these sources 
FWHM MgllA2800 is equal to or greater than FWHM(H/3) and shows a sig- 
nificant blueshift relative to H/3. We interpret the MgllA2800 blueshift as the 
signature of a radiation-driven wind or outflow in the highest accreting quasars. 
In this interpretation the MgllA2800 line width - affected by blueshifted emission 
- is unsuitable for virial mass estimation in « 10% of quasars. 

Subject headings: quasars: general — quasars: emission lines — quasars: indi- 
vidual (SDSS J150813.02+484710.6) 



1. Introduction 



Estimation of black hole mass (Mbh) and Eddington ratio (L/L-Edd) for quasars is of 
great interest both to researchers working on models of the broad line region (BLR) structure 
and to cosmologists. Therefore we need a large number of accurate as possible estimates 
over the widest possible range of redshift, source luminosity and line/ continuum properties. 
FWHM(H/3) is the principal virial estimator at low redshifts. While H/3 can be followed into 
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the infrared (out to z ~ 3.7) at least one additional line is needed to provide complementary 
Mbh estimates over the full redshift range where quasars are observed. MgllA2800 is the 
best candidate since it is a low-ionization line and avoids some of the difficulties associ- 



ated with CivA1549 (e.g. 


Netzer et a 


.1120071 


Trakhtenbrot & Netzer 


2012|; 


Dennev 


l2012h 



The advent of the SDSS database makes possible a direct calibration of MgllA2800 us- 
ing sources where both MgllA2800 and H/3 appear in the same SDSS spectra. Only high s/n 
composite spectra can provide an ideal vehicle for such comparisons. In addition one has to 
consid er that source diversity found within the contex t of a formalism like 4D Eigenvector 1 



f4DEl lSulentic et al 



2000 



2007 



ratio L/Ly, M (e.K. IBoroson fc Green 



1992 



Marziani et al.ll2010 ) is large and likely driven by Eddington 



Marziani et al.l 120011 iBorosonl l2002t iGrupd 12004 



Kruczek et al.ll201ll ; iTang et al.ll2012l ) . The 4DE1 allows to discriminate sources whole emis- 
sion line profiles and whose spectrophotometric properties are strikingly different. A notable 
empirically-motivated boundary at low- and moderate luminosity is set by FWHM(H/3) 
^4000 km s" 1 . Sources narrower than this limit show H/3 profiles that are well fit by a 
Lorentzian funct ion, while broader sour c es show prominent re dward asymmetries in their 
H/3 profiles (e.g. IVeron-Cettv et al.lboOlt iMarziani et aDl2003fi . Sources with FWHM(H/3) 
< 4000 km s _1 include Narrow Line Seyfert Is (NLSyls) by definition and are characterized 
by a significant high-ionization outflow, rev ealed by a CiyA1549 blueshift with respect to rest 
frame or to broad low-ionization lines (e.g. 



Gaskell 


1982; 


Tvtler k Fan 


1992: 



19961 ; ISulentic et al.ll2007t IRichards et al.ll201ll). The outflow has been ascribed t o a radiative 



or m agnetically driven wind (e.g. iMurray fc Chiang! Il997l ; iBottorff et al.lll997l ; iProga et al. 
20001 ) . Therefore it is im portant to at least di stinguish between sources that appear to be 



wind- or disk-dominated (IRichards et al.ll201ll ) appl ying the limit at FW HM( H/3) = 4000 
km s -1 that separates Population A and B sources ( ISulentic et al.ll2000l . c.f. ICollin et al. 



2006|). 



A finer subdivision is still needed even with the restriction to Population A. Pop. A 
sources span a relatively large range in Eddington ratio, AlogL/LEdd~0.5, that likely in- 
volves not only the highest Eddington radiators. We apply the spectral classification of 
Sulentic et all J2002[ ) that divides the plane FWHM(H/3)- J R FeII into bins of Ai? FeII = 0.5 and 
AFWHM(H/3) = 4000 km s _1 , where -RfcII is computed as the equivalent width (or inten- 
sity) ratio of the FellA4570 blend and broad H/3 (see their Fig. 1). Extreme Pop. A (A3 
and A4) sources with -RfcII ~ 1 are ch aracterized by the strongest high-ionizatio n outflow, 
with the largest CivA1549 blueshifts (IMarziani et al.l l2006t ISulentic et al.l 120071 ). and are 
believed to be the highest Eddington ratio sources. A large low- z sample covering the H/3 
and MgllA2800 emission lines is defined (§2]). We point out intriguing changes of the H/3 and 
MgllA2800 line profile occurring in bins A3 and A4 (£J3]) and discuss first order considerations 
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about the physics involved (SJ). 



2. Sample Selection and Construction of Composite Spectra 



We searched SDSS-DR7 for sources catalogued as Type 1 AGN (quasars) in the redshift 
range 0.4 - 0.75 and wit h magnitudes brighter than g »s 18.5 in the g, r or % bands, as well 
as the IZhou et al.l ( 120061 ) catalog. The resultant sample consisted of 716 quasars reduced to 
680 (all of Pop. B: 369 sources; Al: 97, A2: 156, A3: 43, A4:15) by discarding very noisy 
spectra and some sources with unusually red colors. Broad Absorption Line (BAL) QSOs 
were excluded from the sample. 

The rest frame was set by measuring the wavelengths of three of the most prominent 
narrow lines ([Oil] A3727, H/3, and [Qui] A5007 ) when they were detected. Residual sys- 
tematic wavelength shifts ( jHewett &: Wild! 120101 ) in addition to the SDSS-provided redshift 
values were computed taking an average of the three lines in each source spectrum, clipping 
individual meas urements in ca ses of disagreement because of poor data or intrinsic blueshift 
of [Olll]A5007 JHu et alJbooi l. We used IRAF splot to estimate F WHM(H/3) in order t o 
better separate sources into spectral bins following the prescription in lSulentic et al.l ( 120021 ). 
Assignments for all bins were made by visual inspection of each spectrum and estimation 
of .RfcII through ngaussfit. Median composites were constructed respectively for H/3 and 
MgIlA2800 after redshift correction and continuum normalization at 5050 A and 3050 A. The 
rest frame radial velocity of the composites (defined by the average of peak radial velocity 
of narrow [Oll]A3727, E5, H 7 , H/3, and [Om] AA4959,5007) was found < 10 km s" 1 in Al 
and A2, and < 20 km s" 1 in bin A3 and A4, with a rms value always less than 50 km s -1 . 
A line can be considered unshifted with respect to rest frame if \Av T \ < 100 km s" 1 . The 
relative uncertainty of v T measurements for H/3 and MgllA2800 peak velocities on the com- 
posite spectra is somewhat less and has been estimated by propagating 3 sources of error 
on: (1) zero-point; [Oll]A3727 and H/3 narrow component peak wavelengths measured in the 
four composites agree within rms ±10 km s _1 ; (2) wavelength calibration, by measuring the 
dispersion in wavelength measurements of strong sky lines, typically ±15 km s . Systematic 
shifts are consistent with ( £j 25 km s" 1 ) for Hgl 14359 and OH 18401 lines whose wave- 
lengths correspond roughly to the wavelengths of the redshifted MgllA2800 doublet and H/3 
line, respectively; (3) peak line position determined by the multicomponent fit, as provided 
by the fitting program SPECFIT 



I nterpretation of the H/ 3 spect ral range closely follows previous work ( IBoroson fc Green 

19921 : iMarziani et al.l l2003al 120091 ). Along with continuum and Fell emission fit over the 



spectral range 4430 - 5510 A, we will include possible contribution of He I lines at 4471 
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ad 5016 A (that appear to be significant only for bin Al), HellA4686 and a contribution 



been studied by several authors since the mid-1980s ( 


Wampler 


1985 


Graham et al. 


1996; 


Laor et al. 


1997; 


Vestergaard & Wilkes 


2001). 



MgllA2800 are mainly due to Fell. We define a range for specfit analysis (2600-3050 A) 
that is a compromise between proximity to the line and the necessity of having sufficient 
continuum coverage to properly map the broad Fell blends. We also include other known 
l ines e .g. semi-forbidden Al II] 2669.95 and O III 2672.04. in the fits. IBriihweiler fe Verner 
( 120081 ) provide FElluv emission templates co mputed from cloudy simu lations and using 
an 830 level model of the Fe + ion. Use of the iBriihweiler fe Vernerl (120081 ) FElluv template 
results in a systematic residual near 2950 A (Figs. [Hand E]). The excess flux is probably due 
to the blend of He H2945 and Fel emission from a cluster of lines produced by transitions 
from the terms z 5 F and a 5 - D° to the ground state (a 5 D). Fel emission has been predicted 
by p hotoionization models (ISigut et al.l 120041) and was suggested by previous observations 



[e.g.. LK wan et al.1119951 ; IGraham et al.lll996l ). The flux deficit is larger than the Fel predicted 
by photoionization models, at least by a factor o f several with res pect to a low-ionization, 
3 times solar metallicity case (model u20h11 of ISigut et al.ll2004l ). However, several lines 
associated to Fel multiplets 1 and 9 and opt 30 have been convincingly identified in a strong 



Fell emitter ( IGraham et al.lll996f). and the emission of Fe I with respect to Fe II might 



increase with metallicity ( ISigut et al.l 120041 ) . The issue of Fe I emission deserves further 
investigation not last because Fe I intensity is strongly dependent on the assumed Fell 
model at 2900 - 3000 A. For the sake of present paper we checked that the peak shift of the 
MgllA2800 is basically unaffected even by strong changes in the assumed Fe I strength. 



Results 



3.1. Broad Line Profile Analysis 



Composite spectra were analyzed using specfit w ith x 2 minimization techniques ap- 
propriate for non-linear multi-component fits (jKrisslll994f ). The procedure allows for simulta- 
neous continuum, Fell and narrow line fittin g . Two FElluv emission templates were applied: 
the theoretica l one b y IBriihweiler fe Vernerl ( 120081 ). and an empirical template produced by 
Tsuzuki et al.l ( 120061 ). Use of two independent templates was justified by possible effects that 
FElluv subtraction might have on the measurements of line shifts. Continuum subtracted 
composite spectra are shown in Fig. [T]for bins Al and A2, and in Fig. [2] and Fig. [3] for bins 
A3 and A4. The MgllA2800 doublet was first modeled as two Lorentzian-like functions of 
same width and relative intensity ratio 1.25:1. This assumption is justified by the value of 
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the p rototypical A3 source I Zwl, and by the physical conditions within the BLR f lLaor et al. 
19971 ). It is possible that the lines become fully thermalized at the extreme optical depth of 
the low- ionization line (LIL) BLR, justifying the assumption of a 1:1 ratio. We have carried 
out several fits also for the 1:1 case but the results on line shifts and widths are very close 
to the 1.25:1 case and do not affect any of the conclusions discussed below. 

Table [1] reports the results of the our multi-component analysis: intensity, shift ± 
uncertainty at 2cr confidence level (estimated as described at the end of §2]), and FWHM 
of the two line components used to model line profiles: the broad component (BC) and, 
when appropriate, a blueshifted component (blue). The radial velocity was measured with 
reference to the vacuum wavelength of the 2 Ps, — > 2 Si component, 2796.35 A. Uncertainties 

2 2 

in BC centroid shifts are significantly larger toward the line profile base than at peak, and 
are estimated to be ± 400 km s _1 for the centroid at 1/4 peak intensity. Columns BC 
and blue Int. report H/3 and MgllA2800 line intensity normalized by continuum flux at 5050 
A. The normalization at 5050 A ensures that the H/3 value roughly corresponds to the line 
equivalent width, and that meaningful intensity ratios MgllA2800/H/3 can be computed from 
the values reported in the Table. Note that MgllA2800 intensity values are for the doublet, 
while reported FWHMs are for an individual component. The formal uncertainty in FWHM 
measurements of H/3 and MgllA2800 BC is, in absence of systematic effects, around 2% 
for spectral types Al and A2, and ~ 5% in all other cases. Shifts and width of the blue 
components are subject to considerable larger uncertainties since they are close to the much 
stronger BC. In addition their values depend on the line profile assumed for the fit. Formal 
uncertainty (i.e., without considering the possibility of different profile shapes) derived on 
the FWHM of the blue components is ^ 10% in all cases. 

Fits to H/3 and MgllA2800 in Al and A2 bins needed only symmetric, unshifted Lorentzian 
line components. The BC accounts for the entire H/3 and MgllA2800 profiles in Al and A2 
sources where the ratio FWHM(MgllA2800) / FWHM(H/3)^ 0.75 - 0.80 (Ta b. E). This 



FWHM ratio also holds for Pop. B sources and, therefore, for 90% of quasars (I Wang et al. 
20091 : iTrakhtenbrot fc Netzerll2012l . Sulentic et al. 2012, in preparation). The simplest inter- 



pretation is that the emissivity weighted distance of the MgllA2800 emitting gas is somewhat 
larger than that for H/3(«i 1.5 following the virial assumption, see also §1]). 

The H/3 profiles in A3 and A4 bins also involve an (almost) unshifted, symmetric com- 
ponent (the broad component BC) with FWHM(H/?bc) ~ 2000 km s" 1 which we assume 
to be the virial broadening estimator. However, fits to the H/3 profile for sources in A3 
and A4 bins require an additio nal blueshifted (Column "blue" in Tab. ED) component in 



order to minimize residuals (e.g. iLeighlyi |2000| ; iLeighly fc Moordl2004t iMarziani et al.l 12010 



Wang et al.ll201ll ). The "blue" component has been modeled first as a symmetric Gaussian, 
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in line with past work. The right panels of Fig. [2] show that a good fit to the MgllA2800 
profiles in A3 and A4 is possible using a shifted symmetric Lorentz function with a profile 
shift of a few hundred km s _1 . 



At a second stance a skewed Gaussian (jAzzalinilll985l ) has been considered for both H/3 
and MgllA2800 (Fig. [3]). In this case, a two-component model is possible also for the spiky 
MgllA2800 profile. We assumed unshifted BC MgllA2800 emission with FWHM(MgllA2800) 
= 0.8 FWHM(H/3), plus an additional blueshifted component described by a skewed Gaus- 
sian as for H/3. The resulting MgllA2800 line decomposition is shown in the right panels of 
Fig. [3j Line parameters are reported in Tab. [I] (A3b and A4b). 



3.2. A systematic MgllA2800 Blueshift 

Bins A3 and A4 (10% of all quasars) behave differently than the wide majority of quasars 
since their FWHM(MgllA2800) > H/3. The ratio FWHM(MgllA2800)/FWHM(H/3) is larger 
than unity in bin A4, with H/3 and MgllA2800 showing the same width in bin A3 (lower panel 
of Fig. [5]) where a peak blueshift is already highly significant. As mentioned, the MgllA2800 
doublet appears blueshifted with respect to rest frame and H/3 in sp ectral types A3 to A 4 



where H/3 shows evidence for a blueshifted component (Fig. [2} see also lMarziani et al.ll2010l ). 
The MgllA2800 blueshift reaches « 20% of the half line width in bin A4. In bin A3 and 
A4 the core of the MgllA2800 profile is narrow enough to appear visually displaced relative 
to the rest frame (Fig. [2]). Even measuring the position of the broad-line core (without 
any correction because of contaminant lines) with splot we obtain consistent values. The 
effect is too large to be ascribed to sources of uncertainty on rest frame, and it is even more 
significant if relative line shifts are considered. As mentioned, we repeated the fits for all 
bins assuming that the doublet ratio is 1.0:1.0. This results in a slightly larger MgllA2800 
blueshift. To further test the reality of the shift, we considered that the maximum doublet 
ratio for expected physical conditions in the BLR is 1.5:1.0. We constructed noiseless mock 
profiles to derive a peak wavelength in case MgllA2800 treated as single line, for doublet 
ratios 1.5,1.25,1 to 1. The effective wavelength of the doublet is 2799.1, 2799.4, 2800.1 for 
the three ratios respectively. Also in the case 1.5:1.0, the peak shift will remain significant. 

The MgllA2800 A3 and A4 fits with a shifted symmetric function probably yield only 
a lower limit to the shift amplitude since, if a two component interpretation is correct, they 
include also unshifted emission line gas. Yet, these MgllA2800 fits are meaningful since 
they provide a robust measurement of a significant blueshift affecting the MgllA2800 line 
profile. Renouncing the symmetric Gaussian approximation for the blueshifted emission 
(in H/3) provides support for profile decomposition into two components (Fig. [3]). With the 
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exception of A4 H/3, both H/3 and MgllA2800 median spectra show blue components with 
strong blueward asymmetry: the blue co mponent profil es vaguely resemble the "trapezoidal" 
shape of the I Zw 1 CivA1549 profile (c.f. iLeighlylboOoh . The lower v r derived for MgllA2800 
blueshifted component with respect to the one of H/3 in bins A3 and A4 is consistent with 
the profile shape difference since H/3 is more affected toward the line base, while MgllA2800 
is affected to closer to the line core. Indeed, in all cases the shift of H/3 is larger than 
the shift of MgllA2800. This result provides an important constrain on the emitting region 
structure O- Considering the strongly skewed blue component line profiles in A3 and A4, 
the estimated shifts reported in Tabfj] could be considered more properly as upper limits 
since symmetric blueshifted Gaussians would yield some emission under the BC with lower 
shift values. 

The most extreme source in our sample involves SDSS J150813. 02+484710. 6 whose 
MgllA2800 and H/3 profiles are shown in Fig. HI H/3 shows a prominent blue asymmetry 
while MgllA2800 is fully blueshifted with width and shift amplitude similar to the ones 
measured for CivA1549 in extreme Pop. A sources such as I Zw 1. For this source we apply 
to the MgllA2800 profile only the profile model appropriate for CivA1549 of extreme Pop. A 
quasars i.e., an unshifted double Lorentzian + a blueshifted component approximated with 
a skewed Gaussian. While uncommon, other sources like J1508+48 have been found (see 



e.g., Q1258+1404; Barthel et al.lll990h . 



3.3. Major trends 

3.3.1. Spectral types 



The upper panel of Fig. [5] shows radial velocity v T trends of broad component peaks as a 
function of spectral type, while the lower panel shows trends of the ratio FWHM (MgllA2800) 
/ FWHM(H/3). MgllA2800 line shift measurements shown in Fig. 0were carried out in three 
different ways: (1) using SPECFI T with the theo r etical FEII UV template, (2) using SPECFIT 
with the FElluv template from iTsuzuki et al.l (120061 ) and (3) measuring the position of 
the broad-line core without any correction. The three sets of measures yield consistent 
trends. The interesting change in the lower panel of the figure involves a tendency for 
the FWHM ratio to increase fi rst to parity in bin A3 and finally to FWHM(MgllA2800) 
> FWHM(H/3) in bin A4 (cf. iTrakhtenbrot & Netzerl hoiA This trend likely accounts 
fo r the large scatter (and convergence toward parity) of single-source measures in Figure 2 
of IWang et al.l 120091 ). The consistent behavior of both shift and FWHM strengthens our 
confidence that the trends are real. The increase in FWHM(MgllA2800) lends support to 
the hypothesis that an additional blueshifted component is emerging on the blue side of 
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an unshifted BC with FWHM(MgllA2800)^ 0.8 FWHM(H/3). A symmetric blueshifted 
Lorentzian model is helpful for ascertaining the reality of the blueshifts but seems physically 
unrealistic. The blueshifts are more likely associated with a MgllA2800 component due to 
gas moving at larger velocity than the ones inferred from the global shift of the line core 
(Tab. HJ. 



3.3.2. Eddington ratio 



We estimated median Mbh values computing median 5100 A luminosities from fluxes 
of all sources in e ach bin and with F WHM measures of the median composites following 
the prescription of lAssef et al. ( 2011). A bolometr ic correction to the 5100 A luminosities 
was applied following iNemmen fc Brothertonl ( 120101 ) in order to derive L/L-^dd values. The 
semi inter-quartile range (SIQR) of L/L^dd has been estimated using the individual FWHM 
H/3 measurements carried out for spectral bin assignment (§2J). All bins show very similar 
median bolometric luminosities logL « 46.2 [ergs s" 1 ]. Fig. E] shows the peak shift of the 
MgllA2800 broad and blue components as a function of L/L-^dd- Shift values have been 
normalized by the H/5 half-width at half-maximum (HW) in order to provide an indicator 
of dynamical significance to the shifts. The peak MgllA2800 shift is consistent with zero for 
spectral types Al and A2 but begins to appear at log L/L Edd ~ -0.5 at type A2 increasing 
to 0.2/0.3 of half width for A3 and A4. The shift amplitude is much larger when the 
blueshifted component is considered. If 0.8 • HW can be considered as a rough estimator 
of the virial velocity of the MgllA2800 emitting gas then A4 and SDSS J1508+4847 show 
outflows close to escape velocity. The main difference between median A4 and the source 
SDSS J1508+4847 is related to the amount of outflowing gas: in A4 the blue component is 
contributing 1/7 of the total line emission while it exceeds 1/3 in SDSS J1508+4847. 



4. Discussion 



The MgllA2800 profile in median spectra is different from both CivA1549 a nd H/3. Com- 
paris on between the latter two lines yielded inferences about the BLR structure (IMarziani et al 
19961 ). The CivA1549 blueshift in low z sources is thought to be associated with a wind com- 
ponent whose prominence increases with L/L-^dd along the 4D eigenvector 1 (4DE1) sequence. 
The CivA1549 profile in bins A2-A4 can be modeled as a combination of the H/3 BC profile 
plus a fully blueshifted component that acc ounts for most of the flux. This basic scenario 
has been confirmed by recent work (see e.g. iRichards et al.ll2011t IWang et al.ll201lf ). In the 
case of MgllA2800 we see a v T displacement that is also significant with 200 - 300 km s _1 if 
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the shift is measured on the full profile. As pointed out, this is likely a lower limit. If the 
blue component is considered, the shift amplitude is much larger (~ -1000 km s _1 ) but sig- 
nificantly lower in MgllA2800 than in H/3. We are considering median spectra so the profiles 
represent the median behavior of line profiles - in some sense equivalent to a single source 
rms profile. So we can ask: what is the typical relation between the blueshifted CivA1549 
and H/3 on the one hand and the blueshifted MgllA2800 emission on the other? 



4.1. MgllA2800 blueshift: emission from a radiation-driven wind 

The MgllA2800 blueshift is most straightforwardly interpreted as due to outflow motions 
of the line emitting gas with preferential obscuration of the receding part of the flow. The 
fact that large blue shifts are observed when L/ L-^dd is highest indicates a role of radiation 
force in acceleration of the gas. If we consider gravitation and ionizing radiation as the 
only forces (neglecting drag forces and pressure gradients) the radial acceleration can be 
written as a{r) oc crM B n/r 2 [a(r)/(crA r c )L/L Edd — k], where a(r) is the fraction of bolometric 
luminosity absorbed, N c the column density and a the Thompson scattering cross-section, 
and k a constant term. If the first term in square brackets exceeds k, the outflowing velocity 
field of the gas will follow the form v(r) = v t ^Jl — r miD /r, where r min is the launching 
radius of the wind, and v t the terminal velocity (oc y/L/r m i Q ). Outflows driven by line 
and/or ionizing photon pressure can accelerate the line-emitting gas to v t ~ k! (AiL/r)^ ~ 
(M. LJ LvMr Vkp. v \ where t >K Pp i is the Keplerian velocity at the launching radius of the wind 
(e.g. lLaor fc Brandtl 120021 ). In the case of outflows driven by ionizing radiation the force 
multiplier M. is expected to be: M. = ~ 7.5anj,N~^ 3 where a ~0 .5 for Compton-thin 



gas optically thick to the ioni zing; continuum flNetzer fc Marzianill2010l ). Ai can be 3> 1 in 
the case of line driven winds ( Proga et al. 200ol ). 



The largest CivA1549 blue shifts are observed in spectral types A3 and A4 (in A3, 
— 1000 ^ Av T £ — 2000 km s" 1 ) where also the MgllA2800 shifts occur. It therefore seems 
unlikely that there is no connection between CivA1549 and MgllA2800: both lines may be 
emitted as part of the same flow. However, bulk emission are expected to occur at differ- 
ent distances and/or in different physical condition. In the framework of photoionization, 



f Netzer 


1980; 


Korista et al. 


1997) 



mainly beyond the fully ionized zone of geometrical depth h ~ 10 23 £/riH _1 (with U being 
the ionization parameter) where all of the CivA1549 is emitted. The total column density 
needed for substantial MgllA2800 production is not well constrained at very low ionization, 
since the fully ionized zone is already a tiny fraction of the emitting gas slab if iV c ~ 10 23 
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< 



cm and log U 

Since v t is proportional to both N c ~ 1 / 2 and r~V 2 , MgllA2800 emission may occur at 
higher column density and/or larger distance than blueshifted CivA1549. A large r min is 
consistent with the overall symmetry of the MgllA2800 profile base ( §3.2( Figs. HI El [3j) . 
Reverberation mapping indicates that hi gh ionization lines are em itted closer to the central 



continuum than low- ionization lines (e.g. iPeterson fc Wandellll999l ). It is possible to ascribe 



the blueshifted MgllA2800 emission entirely to a larger radial distance if MgllA2800 arises 
~ 2 times more distant than CivA1549. On the other hand, if L/L Edd — > 1, v t ~ VA-Ttw- 
Restricting our considerations to order-of-magnitude estimates due to uncertain shift values, 
a reasonable increase in column density over the standard value A c ~ 10 23 cm" 2 would also 
suffice to reduce the MgllA2800 shifts close to the observed values. 

An assumption from the behavior of the FWHM ratio is that the MgllA2800 flows start 
close to where the bulk of the low-ionization lines (H/3 in the present case) is emitted. The 
following considerations apply also if FWHM(MgllA2800)/FWHM(H/3)^ 0.8 implies «1.5x 
larger emissivity- weighted distance for MgllA2800 than for H/3. The distance from the 
continuum source of the emitt i ng ga s r em can be derived from the continuum luminosity at 



5100 A following iBentz et al.l ( 2009 ): logr em « 17.6 [cm] for spectral type A4. If logriH= 
12 [cm -3 ], the ionization parameter is \ogU ~ —2.8. The MgllA2800 luminosity is then 
L(MgllA2800) = 47rr 2 m / c cr em = 47rr 2 m / c e^-, where f c is the covering factor, a p m the emerging 



line fl ux for unit surface, and e the depth-averaged volume emissivity. cloudy ( IFerland et al. 



19981 ) simulations indicate that a cm has a minimum value ~ 10 7,8 ergs s 1 cm 2 at log ran = 
12 if -2.8 < log U < -2, 10 23 cm" 2 < N c < 10 25 cm~ 2 , 12 < logn H < 13. In these ranges a em 
depends slightly on A c and U. The computed L(MgllA2800) extrapolated to full continuum 
coverage is always larger than the observed L(MgllA2800) (~ 10 43 ' 5 ergs s" 1 for A4). This 
is the case if the geometry is assumed static and open or if a velocity field appropriate 
for a wind (i.e., with photon local escape probability following Sobolev's approximation) is 
considered. The derived f c ^5 0.2 indicates partial covering of the continuum as in a wind or 
in an ensemble of outflowing clouds. 

Resonant line acceleration is expected to contribute to the dynamics of the flow in the 
physical scenario outlined above. Circumstantial evidence in favor of line acceleration is 
provided by the difference between a resonance UV line (MgllA2800) and a non-resonance 
line (H/3). From the purely observational point of view, a line driven outflow would be 
convincingly demonstrated if the ionizing photon flux were found unable to drive the line 
emitting gas to the observed outflow velocity or, in the context of absorption lines, from 



"line locking" (e.g. ICottis et al.ll2010l ). However, the ionizing photon flux appears to be 



sufficient to accelerate the gas to the observed MgllA2800 velocities and to escape velocity 
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in A4 sources: for logL/L E dd~ -0.2, M « 7.5, v t « 1.8v K epi- If ^'Kcpi ~ 0.8HW(H/3), 
f t ~ 1500 km s" 1 . This value exceeds the peak velocity of the blue component in bin A3, 
and is in agreement with the ones measured in A4 and in SDSS J1508+4847. Resonant line 
acceleration might be needed if the gas has a large N c ( £ 10 23 cm -2 ) or if the MgllA2800 
emitting gas is shielded by part of the continuum. 



4.2. Alternate interpretations 



Both MgllA2800 and H/3 line profile widths are probably modified by the viewing angle 



of the outflow/jet axis. Evidence exists that the 



orientation of the 


et axis 


1995 




Rokaki et al. 


2003 





ine width of H/3 ' is aff e cted by line-of-sight 



20031 : ISulentic et al.ll2003l : IZamfir et all 120081 ). An effect on FWHM of a 



factor ^2 is likely between core an d lobe-dominated sources. Orientation effects a re also ex- 
pected for radio-quiet quasars (e.g. lJarvis fc McLurdl2006l : iPunsly fc Zhangll2010l ). Extreme 
and variable soft X-ray emission from some narr ow line Seyfert Is (P op. A) sources has been 
interpreted as a signature of pole-on orientation ( ISulentic et al.ll2000l . and references therein). 
Recent work confirms a dependence on orientation for H/3 in radio- lou d sources and further 
suggests a less-strong dependence for MgllA2800 (IRunnoe et al.ll2012l ). Following this line 
of reasoning the occurrence of blueshifts in bin A3/A4 might involve sources viewed at a 
favorable line-of-sight orientation. It is not clear whether the results of iRunnoe et al.l (120 12[ ) 
on different orientation sensitivity can be extended to Pop. A where radio-loud sources are 
almost absent in bins A2, A3 and A4. Even if the FWHM(H/3) change could be explained on 



the basis of an orientation e 
Al to A4 



Wi 



2004; Baldwin et al. 



s et a. 



2004 : 



ect, several line intensity rat i os change very strongly going from 



1999 



Aoki fc Yos 



Negrete et al. 



2012 



lida 



1999 



Shin et al. 



Sulentic et al. 2000 



Bachev et al. 



20121 ) : -RFeii by definition, but also 



Cm] A1909 /Sim] A1892, AliiiA1860/Siin]A1892, SiivA1397+Oiv]A1402/CivA1549. Emission 
line equivalent widths change as well. For instance, the EW of H/3 BC shows a de crease from 
Al to A4 by a factor ^2 (Tab. [[]), as found previously (e.g., ISulentic et al.ll2000l ). The large 
difference in EW persists also if the flux of the blueshifted component is included. Line 
intensity ratios and line equivalent widths are most likely sensitive to density, ionization 
state and chemical composition of the gas along with ionizing continuum shape (ultimately 
thought to be governed by Eddington ratio). In this respect we note that the decrease of 
H / (MgllA2800) from ^110 to m 70 A is als o consistent with a study showing an anticorre- 
lation between iy(MgllA2800) and L/L-sd& (IDong et al.ll2009l ). It is unclear how orientation 
might drive such changes. We are probably dealing with a restricted range of L / L-fm in each 
spectral bin "convolved" with the effect of orientation f lMarziani et al.l 120011 ) . If a (rare) 
pole-on orientation favors the observation of large shifts then SDSS J1508+4847 might be 



an example of a pole-on source. 



5. Conclusion 



MgllA2800 should be used as a virial estimator with caution in high Eddington ratio 
sources. Under the simplest assumptions virial motion implies an unshifted and symmetric 
profile. A systematic line profile blueshift can be interpreted as the signature of emission from 
radiatively acceleration of gas motion therefore invalidating the virial broadening assumption 
for MgllA2800 in 20% of Pop. A sources (10% of all quasars). Conversely further work 
has shown that the majority of quasars show unshifted MgllA2800 profiles that are more 
symmetric than H/3 (Sulentic et al. 2012 in preparation, and references therein). The width 
of MgllA2800 is probably a suitable virial broadening estimator for those sources. 
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Table 1: Derived Quantities from the EL3 and MgllA2800 Profile Multicomponent Analysis 



Sp. 


T. 








H/3 










Mj 


;llA2800 








BC 






blue 






BC 






blue 




Int. a 


Shift b 


FWHM b 


Int. a 


Shift b 


FWHM b 


Int. c 


Shift b 


FWHM 


b Int. c 


Shift b 


FWHM b 


Al 




96 


20±40 


3180 


0.0 


n.a. 


n.a. 


111 


35±50 


2710 


0.0 


n.a. 


n.a. 


A2 




85 


-20±40 


2900 


0.0 


n.a. 


n.a. 


73 


-70±50 


2320 


0.0 


n.a. 


n.a. 


A3 




15 


45±40 


2190 


12 


-1240±230 


4250 


65 


-150±50 


2240 


0.0 


n.a. 


n.a. 


A3b 




51 


45±40 


2190 


6.5 


-1420±510 


4100 


54 


d 


1750 c 


5.6 


-880±170 


3100 


A4 




28 


70±50 


1980 


16 


-1240±220 


4870 


68 


-265±90 


2650 


0.0 


n.a. 


n.a. 


A4b 




29 


70±50 


1940 


13 


-1530±200 


4460 


46 


d 


1585° 


10 - 


1010±130 


3300 


J150J 


S+48 


32 


-5 d 


2300 


5.0 


-1540±100 


4000 


37 


d 


2190 


20 - 


1490±160 


3500 



Line intensity normalized to the continuum at 5050 A. The value roughly corresponds 
to the rest-frame equivalent width in A. For J1508+48 values are in units of 10~ 15 ergs 
s -1 cm -2 A -1 . 



b In units of km s 1 . 

c Line intensity normalized to the continuum at 5050 A. The value can be used an 

estimate of the MgllA2800/H/3 intensity ratio. For J1508+48 values are in units of 
IQ-15 er g g s -i cm -2 A -1 . 

d Imposed to be consistent with rest frame. 



e FWHM(MgllA2800) = 0.8 FWHM(H/3), for broad component. 
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Fig. 1. — Spectra of H/3 (left panels) and MgllA2800 (right panels) for spectral types Al 
(top) and A2 (bottom). The horizontal scale is rest frame wavelength [A] or radial velocity 
with the origin indicating rest frame (laboratory) wavelength. In the MgllA2800 panels the 
vertical dot-dashed line is drawn at the reference wavelength 2799.4 which corresponds to 
a component ratio 1.25:1.00. The black lines show the original continuum-subtracted spec- 
trum while the dashed magenta line shows the model including all emission line components. 
The thick black lines show the broad component and thin black lines the individual com- 
ponents of the MgllA2800 doublet. The green lines trace Fe ll opt and FElluv emission and 
the gold- brown lines various contributions associated with the narrow-line region (H/Snc, 
[Oiii]AA4959,5007). In the H/3 panels significant HellA4686 is revealed (thick blue line). In 
the Al H/3 panel He H4924 and He 115016 (black lines) almost overlap with the m42 Fe ll opt 
lines. The Fel + He 112945 emission (modeled as the sum of two Gaussians) is also traced 
by a brown line and is visible toward the right end of the MgllA2800 panel at ~ 2950A. 
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Fig. 2. — Spectra of H/3 (left panels) and MgllA2800 (right panels) for spectral types A3 (top) 
and A4 (bottom). The right panels show a model that assumes a shifted BC for MgllA2800 
and an almost mishitted BC + a blueshifted Gaussian component (thick blue line) for H/3. 
Meaning of all other symbols is the same as Figure HJ 
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Fig. 3. — Spectra of H/3 (left panels) and MgllA2800 (right panels) for spectral types A3 
(top) and A4 (bottom), as in the previous figure. The panels show models of both H/3 (left) 
and MgllA2800 that assume an (almost) unshifted BC and a blueshifted skewed Gaussian 
component. Meaning of symbols is the same as Fig. |2] except for the individual MgllA2800 
components that are shown by thick lines in the right panels. Refence dot-dashed lines are 
drawn at the vacuum laboratory wavelength of the two MgllA2800 components. 
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Fig. 4. — Spectra of H/3 (left panels) and MgllA2800 (right panels) for quasar SDSS 
J150813. 02+484710. 6. Vertical scale is specific flux in units of 10" 15 ergs s" 1 cm" 2 A" 1 . 
The meaning of all other symbols is the same as Figs. [U and [3j 
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Fig. 5. — Trends in Pop. A as a function of spectral type. Filled circles: MgllA2800; open 
circles: H/3. In the upper panel the BC shifts are relative to rest frame. Error bars are at 2a 
confidence level, and are shown only for one measure of H/3 and one of MgllA2800 for clarity. 
The squares indicate an additional measurement of the core centroid of the line obtained after 
removal of NC emission. The re d circles refer to measurements obtained applying the Fell 
template of iTsuzuki et al.l ( 120061 ) in the region of MgllA2800. The dot dashed lines at ±100 
km s _1 in the top panel defines a zone where peak measurements are consistent with zero 
shift (at 2cr confidence level). In the lower panel the ratio FWHM(MgllA2800) (individual 
component) over FWHM(H/3bc) is plotted. Error bars are again at la confidence level. 
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Fig. 6.— The ratio between the peak MgllA2800 shift and H/3 half-width half maxi- 
mum (HWHM) as a function of Eddington ratio. Data points are shown for an unresolved 
MgllA2800 fit (Al, A2, A3, A4), as well as for the blueshifted component in two component 
fits (A3b and A4b). Error bars are at 2a confidence level. SIQRs in L/L^dd values are ^0.1 
for bins Al,A3,A4, and ~ 0.16 in bin A2. Eddington ratio ranges corresponding to ±SIQR 
are indicated in the lower part of the diagram for all four spectral types. The blue data 
point refers to a two component fit for source SDSS J150813. 02+484710. 6. 
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